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^ 1 Abstract 

We examine the concentration invariance of cyclic species for the irreversible polymerization. The simulation re- 
£f*^ | suits show that the invariance holds in good approximation in the irreversible process also. The validity of the 
simulation equations employed was confirmed through comparison with the recent experimental observations by 
the Nagoya university group carried out independently of the present work. 
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O ! 1 Introduction 

1*3 i 

It has been well established that the relative frequency of ring formation to intermolecular reaction increases 
with dilution. It was pointed out, on one hand, that the ring concentration is independent of the monomer 
concentration. This feature has been proven rigorously for the reversible difunctional system [1]. In this paper we 
examine the irreversible polymerization [2] to prove that the concentration invariance of cyclic species is satisfied for 
the irreversible case also. The soundness of the result is confirmed through comparison with the recent experimental 
observation of the telechelic polystyrene by the Nagoya university group [3] . 
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ID ! 2 Theoretical 

Let Mo be the initial monomer number and V the system volume. The birth-death equation for the AA-typc 
difunctional polymerization can be expressed in the form: 
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where N x and Nj^ denote the numbers of chain x-mers and cyclic x-mers respectively (x = 1,2, ■■ ■ being the 
number of repeating units), u the number of unit reactions (u = 0, 1,2, •••), and tp x the relative cyclization 
frequency [1]. For the bifunctional system u is equal to the number of bonds. 

The quantity ip x is defined by & /v where & is the probability that one end on a linear x-mer chain enters the 
small volume v with the radius i s around the another end. For the Gaussian chain <p x has the form: 

<p x = {d/2n d 'H d s N A ) t'^ti-^dt, (3) 



where d is the space dimension, Na the Avogadro number and v x the quantity defined by 



(ri)/i a 2 , (4) 



where (r^) is the end-to-end distance. In poor solvents, or in the limit of infinite concentration (C — ► oo), the 
excluded volume effect is expected to vanish and then 

v x = (rl)e/tt=C F Z e x (4') 



with the subscript signifying the Flory unperturbed state, Cp the Flory characteristic ratio, and £ e the effective 
bond number within a repeating unit defined by £ e = Yli H- Note that Cf is not a constant, but an increasing 
function of x and approaches a constant as x — > oo. 

Contrary to the ideal system without rings, there is no explicit closed solution [1] for the above equations, so 
N x or Nr x can not be expressed uniquely in terms of u. Fortunately, eqs. (1) and (2) can be solved numerically [2], 
Chemical reactions are by nature discrete, because the number of events is enumerable. Thus these equations are 
essentially recursive. 

To solve eqs. (1) and (2), it is convenient to sum up the above equations from u = to p (an arbitrary number 
of events) to obtain 
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with the boundary conditions: 

N x {0) = \ ° t X * ] N Rt (0) = for all x's. 
Mo if x = 1, 



Table 1: Model polymerization of the telechelic polystyrene in THF 



parameter 


symbol 


value 


molecular weight 


M w 


4.6 x fO 4 


characteristic ratio 


C F 


10 


expansion factor 


a 


1.39 


effective bond number 


£e 


867 


standard bond length (A) 


is 


1.55 


monomer concentration (mol/l) 


C 


2 x icr 4 2i 



Let us perform the calculation of the above equations modeling the 
polymerization of the telechelic polystyrene (t-PS: Mw = 4.6 x 10 4 , x = 
440) in tetrahydrofuran (THF) as a good solvent with the help of the 
parameters shown in Table 1 [3]. We assume the Gaussian behavior of 
the polystyrene chain having the expansion factor a = 1.39, i.e., the 
Gaussian chain having the end-to-end distance (r^) = a 2 (r^)g). 

The simulation was performed to calculate the weight average dis- 
tribution for D = 0.84 (Fig. 1) and the limiting case of D — > 1 (Fig. 
2), where D is the extent of the advancement of reaction defined by 
D = u/M . 

In Fig. 1 solid lines represent the theoretical curves based on eqs. (5) 
and (6) (D = 0.84 is assumed) for cyclic and acyclic species, respectively, 
and diamonds ({>) and crosses (x) are corresponding experimental obser- 
vations by Kushida and coworkers [3] under the same reaction condition. 
Agreement between the theory and the experiments is amazingly good. 
This convinces us the soundness of the kinetic formulas (5) and (6). 

From Fig. 2, one can see that the concentration of rings ([r]) has the 
T-line curve which was already observed in the equilibrium polymerization 
(see Remark 1); i.e., there is a plateau zone in high concentration where 
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Fig. 1: Experimental weight fractions 
of cyclic (<C>) and acyclic species (x) 
in the polymerization of the telechelic 
polystyrene monomer (MW = 46000) (Re- 
sults from Kushida and coworkers [3]). 
Solid lines (— ): theoretical lines at D — 
0.84 by eqs. (5) and (6). 



t 1 corresponds to = 0.1w/v%. 



[r] is approximately constant over broad monomer concentration (Co) range, whereas in the low concentration of 
Co < 10 -3 mol/l, [r] decreases strongly with increasing dilution along with the sudden rise of the weight fraction 
wr of rings. The upper bound of Co (= 0.02) in Fig. 2 corresponds to the non-solvent system. The strong rise of 
wr in the low concentration indicates that the ring formation begins to overwhelm the intcrmolccular reaction, in 
which the production of large molecules is no longer possible because the molecular growth occurs only through 
intcrmolccular linking. 

As shown earlier, eqs. (5) and (6) has the limiting solution of Cq — * oo: 



As D — ► 1, eq. (7) leads to 
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Fig. 2: Total concentration of cyclic 
species and the weight fraction as against 
initial monomer (PSt) concentration at 
D = 0.99. [F] and wr were calculated ac- 
cording to eqs. (5) and (6); dotted line 
(•••): eq. (8). 



The dotted line in Fig. 2 illustrates this eq. (8). One can see that there is 
slight numerical difference in the plateau zone between [r] of the present 
simulation (solid line) and that of eq. (8) (dotted line). There is also 
observed substantial difference in quality of the concentration invariance 
of [r] between the irreversible case (Fig. 2) and the reversible case (Fig. 
3). The differences is (probably) due to the small sample space (M = 
100) employed in this simulation. The essential feature is, however, quite 
alike, namely, [r] is almost constant over broad concentration range. Thus 
one can apply the invariance principle to the irreversible polymeization 
process without introducing an appreciable error to our final results [1] . 

An important question is whether one can extend the principle of the concentration invariance of cyclic species to 
a more general case of / > 3 [1]. At present a rigorous proof does not appear available. From the chemical dynamics 
point of view, however, the two systems (linear and branching) have exactly the same feature that the molecular 
growth occurs through the competition between the intermolecular reaction and the cyclization of linear molecules 
(see Fig. 3). In this sense, no difference exists between the linear system and the branching system: one may regard 
every reaction in a branched system as a reaction of a chain molecule (cyclization) or 
between chain molecules (intermolecular reaction) embedded on clusters. In fact we 
can describe the ring- forming process of the branching system by the formula precisely 
identical with that of the linear system: 
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For a large Co, it follows that ^2 ( v R^/ v l) "C 1, and the above equation yields the 
known limiting solution: for instance, Vx [(/ — 1)D] X /2x for the R-Ay model [1]. 
The solution is independent of Co, which convinces us strongly the principle of the 
concentration invariance of the cyclic production for the branching process. 
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Fig. 3: Representation of a 
chain on a cluster. 



3 Possibility of Experimental Verification of the Concentration Invari- 
ance 

Let Q denote the total number of molecules in a given branching system and Qq the corresponding quanatity in 
the ideal system with no rings. Following the known relation, we write as: 

ft - n = r, (10) 

while 



X n = M /Q 
n = M (l - \fD) 



(11) 

(12) 



where X n is the number average molecular weight. Eq. (10) is a universal relation valid for all bond forming 
processes within the interval < D < D c . Substituting eqs. (11) and (12) into eq. (10) and using the relation 
C = Mq/V and [r] = T/V, one has the equality: 

i r i= c °(^- i+ H' (i3 » 

The conclusion of eq. (13) is important. It states that if X n can be measured experimentally as a function of D 
and Co, one can estimate [r] as a function of Co- This provides a possibility that one can test experimentally the 
assumption of the concentration invariance of [r] for general branching processes of / > 3. This possibility was 
first suggested by Faliagas [4]. To apply eq. (13), it can be conveniently recast in the form: 

Xn = [ryco + i-ifD- (14) 

From Fig. 2, the term (l — hfD) is of the order « 0.01, while the quantity [r]/Co is of the order w 10~ 3 . Hence 
the effect of the ring formation on X n is experimentally measurable. 

4 Conclusion 

The present simulation showed that the concentration invariance is well satisfied for the irreversible difunctional 
polymerization. There exist some reasons to consider that this principle applies to the branching process also; for 
instance, (i) the limiting solution is independent of Co; (ii) no inconsistency has been found up to present for the 
theory of gelation constructed on the basis of this invariance principle [1] . There is a possibility to prove directly 
the concentration invariance of cyclic species through the relationship of eq. (14). 

Appendix 

It will be useful to compare the result in the text with the equilibrium case [1] . Assume an imaginary case in which 
the above-mentioned PSt undergoes the equilibrium polymerization. The corresponding equilibrium solution is of 
the form: 

\1-wrJ 

x—Q v ' 

so that 

C = ±±^(^-Y/2. (16) 
wr e -~L \1—Wr ' 

x—ij 

We show below a special case of D — ► 1, 
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Fig. 4: Total concentration [F] and weight fraction wr of cyclic species in equilibrium polymerization 
as against initial monomer (PSt) concentration Co. Fig. 4-(B) is a magnification of Fig. 4-(A). 
The arrows (j) in Fig. 4-B indicate the critical concentration, Cq, below which the concentration 
invariance of cyclic species breaks down. 



As one can see, there is a clear-cut critical point, Cg, below which the concentration invariane of [r] breaks down. 
From eq. (16), it follows that Cq = Y^Lo vW^- For Co < Cq, the weight fraction of rings, wr, is held unity, while 
[r] decreases with dilution; first slowly, then rapidly, and finally as [r] = Co, showing that the monomcric ring 
occurs more abundantly with dilution and exclusively in the limit of Co — » 0. 
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